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INTRODUCTION 


A computer code developed previously [1] has been used to simulate lateral compositional 
distributions for several II- VI semiconducting systems. The code is based on analytical results [2] 
obtained by following Coriell et al. formalism [3], The system under study is an azimuthal 
symmetric cylindrical system with a curved melt-solid interface shape during an unidirectional 
solidification process. It is assumed that the system is under a steady state diffusion limited 
growth condition with a curved melt-solid interface and equation for the diffusion limited growth 
is 


V^C (r, z) + p|-C (r, z) = 0 . 
d z 


( 1 ) 


In terms of cylindrical coordinates explicitly, the equation is expressed as 


2 2 

A* (r, z) + C (r, z) + ArC (r, z) + pf (C (r, z)) = 0, 


dr 


r dr 


dz 


dz 


( 2 ) 


where the dimensionless solute concentration, C, has been scaled by the bulk concentration Cq, r 
and z are dimensionless axial and radial coordinates which have been scaled by the radius of the 
ampoule, R, and p=VR/D, in which V is the growth velocity, R is the sample radius and D the 
diffusion constant of the solute in the liquid. 

The boundary conditions to be satisfied are 

C (r, z — » oo ) = 1 , (3) 

■^C (r, z) = 0 at r=l and r=0, (4) 

dr 


P(k-l)Cj(r, 2 ) =^ C I ( r ' ! )-|c [ (r, 2 )|w(r). (5) 

where k is the segregation constant, W(r) the melt-solid interface and I denotes at the melt-solid 
interface location. 

The melt-solid interface shape can be expressed as a series of Bessel’s functions. 
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00 


W(r) =8+ £ 6(n) J Q (u n r) . 
n = 1 


( 6 ) 


where 8 is the planar interface location and 8(n)’s are the coefficients for the Bessel’s functions 
J 0 (u n r) and u n ’s are the zeros of Jj(a). 

The diffusion equation, equation (2), can be solved by the separation of variables method. 
In the limit that dw/dr is small compared to 2k+p(n)-l, where p(n) denote I (2u_)~ and have 



(2u 

+ 

n 

i 

P J 


large values for small P, the solution for solute concentration in the melt at coordinate (r, z) is 

c ( „) . 1 + LzV^- 6 > + y 

( * } k 6 + 2- 2k + p (n) — 1 C ,(?) 

n = 1 

The solute concentration in the melt at the melt-solid interface is 


00 


C,(r,z) = 1 + ' J‘ e -P(W-8) + £ 


n = 1 


2p(l -k)S(n) 
2k + p (n) — 1 


( 8 ) 


The solute concentrations of the solid at the interface, Csi(r, z)=kCi(r, z), is 


00 


C (r,z) = k+ (1 -k)e -p(W_5) + V 2pk(1 k)5(n) 

SI ^ 2k + p(n) - 1 

n = 1 


( 9 ) 


For small P limit and P(W-8) is small, we may expand the exponent term and obtain the 
solute concentration in the solid at the melt-solid interface as 


00 


C gI (r,z) = l-p(l-k) £ 


8 ( n ) Jq ( u n r ) 


n = 1 1 + 


2k 

p(n) - 1 


( 10 ) 
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If 2k/[p(n)-l] « 1, then equation (10) reduces to a simpler expression 


C SI (r,z) = l-p(l-k) [W (r) -6], (11) 

i.e., the deviation of C§] from unity is proportional to the deviation of interface shape from 
planarity. The proportional constant is the product of P and (1-k). 

NUMERICAL CALCULATIONS AND RESULTS 

We use equations (9), (10), (11) to calculate compositional distribution at the melt-solid 
interface for a diffusion limited growth system with a curved melt-solid interface shape for small 
P values. The results in Figure 1 show that equations (9) and (10) give very close answers while 
equation (11) over estimates the results by a few percent. Calculated lateral compositional 
distributions by using equation (10) for four different melt-solid interface shapes are shown in 
Figure 2. Scaling these compositional distributions with their interface deflections yields curves 
with similar shapes as shown in Figure 3. By adjusting the scaled curves they all fall on one curve, 
indicating that the lateral compositional distribution is proportional to melt-solid interface shape 
for same P and k (see Figure 4). The results agree with the simple expression equation (11) 
because both the conditions, P is small and 2k/[p(n)-l] « 1, are satisfied in this study. The 
experimental conditions for various systems under consideration are listed in Table 1 which have 
P equal to 0.1455, 0.16 or 0.2667 small enough such that the analytical and approximate 
expression obtained above apply. In the MZT USML-1 experiments [4], the interface shapes were 
not radially symmetrical, therefore, the deflections of the interfaces for the right and left handed 
sides were measured. Table 2 is a comparison of the calculated and experimental radial 
segregation on MZT USML-1 experiment. The results agree reasonably well. Figure 5 shows the 
simulated and experimental data for MCT D4 crystal grown under axial 50 kG magnetic field [5], 
The experimental data agree quite well with the theoretical results indicating that the convection 
in the melt during the growth has been reduced by the magnetic field. Figure 6 gives the simulated 
and experimental data for MCT D6 crystal grown without magnetic field [5], Due to larger 
convective fluid flow in the melt during crystal growth in D6 experiment compared to that of D4 
experiment, D6 crystal has larger deflection of melt-solid interface shape. The simulated results 
show that the composition distribution scales with melt-solid interface shape, but with larger 
discrepancies between the two curves at the edges. This is due to the contribution from the 
convective fluid flow in the melt. The simulated data and the experimental data for MZT B16-4 
crystal are shown in Figure 7. The disagreement between theory and experiment in MCT D6 and 
MZT B16-4 experiments are also believed to be due to contributions from convective fluid flow 
[6], Figure 8 shows the simulated and experimental data for MCT USMP-2 space flight crystal [7] 
grown in space. The two curves have significant mismatches at the edges. Which was believed to 
be contributed by the residual transverse microgravity. 
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Figure l.Cl,C2 and C22 are calculated 
lateral composition distributions 
for an MCT sample with k= 4, 
(3=0. 1455, and deflection= 0.8125 
by using equation (9), (10) and 
(11) respectively 



0 0.5 1 


Figure 2. Composition distributions for 

various deflections. The radius of 
the sample is 0.4 cm. The deviation 
of the interfaces are 0.092mm, 

0.136 mm, 0.368mm & 0.544 mm. 



Figure 3 Deflection scaled composition Figure 4. Adjusted deflection scaled composition 

distributions. The radius of the sample is distributions. The radius of the sample is 

0.4 cm. The deviation of the interfaces 0.4 cm. The deviation of the interfaces 

are 0.092mm, 0. 136 mm, 0.368mm & are 0.092mm, 0.136 mm, 0.368mm & 

0. 544 mm. 0 .544 mm.The results show the lateral 

composition deviation scale with the 
interface deviation. 
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Figure 5. Calculated and measured data of MCT D4 
crystal with 50 KG magnetic field 


Figure 6. Calculated and measured data of MCT 
D6 crystal without magnetic field 



Figure 7. Calculated and measured data of MZT 
B6*4 ground based crystal. Deflections 
are 1.086 & 1.091 on left and right hand 
sides of the crystal. 



*(mm> 

Figure 8. MCT-27Q, Right Sample Composition 
across x-Diameter 
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Table 1: Ampoule Data 


Ampoule 

I.D. 

C 0 

V (cm/s) 

fl 



P=VR/D 

Deflection 

A 

MCT D4 

0.2 

2x10-5 

m 

5.5x10-5 

0.4 

0.1455 

0.275 

MCT D6 

0.2 

2x10-5 


5.5x10-5 

0.4 

0.1455 

0.8125 

MZT 

USML-1 

0.16 

4x10-6 

m 

6x10-6 

0.4 

0.2667 

R-0.5 

L-0.225 

MZT 

B16-4 

0.16 

4x10-6 

m 

6x10-6 

0.4 

0.2667 

R- 1.091 
L- 1.086 

MCT 

USMP-2 

0.2 

2.2x10-5 

m 

5.5x10-5 

0.4 

0.16 

0.5 


Table 2: Measured and simulated lateral compositional deviations of MZT USML-1 space 

flight experimental data 


Measured 

Calculated 

ACj/qpO.40 

ACi/C 0 =0.375 

AC2/Cq=0.18 

AC 2 /C 0 =0.153 
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